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\b" ' ABSTRACT 

o 
o , 

Theoretical studies predict that Trojans are likely a frequent byproduct of planet 
formation and evolution. We present a novel method of detecting Trojan companions 
to transiting extrasolar planets which involves comparing the midtime of eclipse with 
the time of the stellar reflex velocity null. We demonstrate that this method offers the 
potential to detect terrestrial-mass Trojans using existing ground-based observatories. 
This method rules out Trojan companions to HD 209458b and HD 149026b more massive 
than ~ 13 M ffi and ~ 25 at a 99.9% confidence level. Such a Trojan would be 
dynamically stable, would not yet have been detected by photometric or spectroscopic 
^vq \ monitoring, and would be unrecognizable from radial velocity observations alone. We 

outline the future prospects for this method, and show that the detection of a "Hot 
Trojan" of any mass would place a significant constraint on theories of orbital migration. 

o " 

■ Subject headings: techniques: photometric, radial velocities — planetary systems: for- 

i 1 mation — celestial mechanics 

o 

to 

1. Introduction 

^ Stable Trojan companions to extrasolar planets may be common. In our solar system, Mars, 

Jupiter, and Neptune share their orbit with asteroids orbiting near the stable (L4/L5) Lagrange 
points that lead/trail the planet by ~ 60°. Orbits near the L4/L5 points of the terrestrial planets, 
Saturn and Uranus are less stable due to perturbations from the other planets (Nesvorny & Dones 
2002). Saturn's satellites also include small moons orbiting about the L4/L5 points of Tethys 
and Dione. While the mass ratios of the Trojan systems in our solar system are rather extreme 
(< 7 x 10~ 9 ), extrasolar planets may have more massive Trojans. Theorists have already outlined 
several mechanisms to form Trojans with mass ratios as large as unity. For example, Laughlin & 
Chambers (2002) present hydrodynamic simulations of a protoplanetary disk where disk material 
lingers near the L4 and L5 points of a planet (near the gap-opening threshold). The resulting 
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vortex could trap particles and lead to the accretion of a Trojan in situ (Chiang & Lithwick 2005). 
If disk torques caused the planet to gradually migrate inwards, then the Trojan would migrate 
with the planet. Unlike resonant migration in the 2:1 mean motion resonance, the eccentricity or 
libration amplitude of the Trojan would not be excited by the migration (Laughlin & Chambers 
2002). Alternatively, a body could be captured into an orbit about the L4/L5 point after a violent 
event, as has been suggested for the formation of Jupiter and Neptune Trojans in our solar system 
(Morbidelli et al. 2005). Capture into a Trojan orbit could also occur due to rapid mass growth 
of the planet or a collision of two objects near L4/L5 (Chiang & Lithwick 2005 and reference 
therein). Another possibility is that convergent migration could trap multiple protoplanets into a 
1:1 mean motion resonance (Thommes 2005, Cresswell & Nelson 2006). In each of these scenarios, 
the captured bodies could initially have a large libration amplitude or reside on horseshoe-type 
orbits. However, if the capture occurred before or during the planet's inward migration, then 
interactions with either a gaseous or planetesimal disk would damp the libration amplitude. This 
mechanism is even capable of causing objects initially on horseshoe orbits to evolve into tadpole 
orbits and then small amplitude libration near the L4/L5 fixed point. Such behavior has been 
found in numerical simulations of multiple planet systems interacting with either a gaseous or 
planetesimal disk (Cresswell Sz Nelson 2006; Ford & Chiang 2006). 

Trojans of Jupiter and Neptune have provided clues about our solar system's history (Michtchenko, 
Beauge & Roig 2001; Kortenkamp, Malhotra & Michtchenko 2004; Chiang k Lithwick 2005; Mor- 
bidelli et al. 2005). Similarly, the detection of extrasolar Trojans would be useful for constraining 
theories of planet formation. While all the above mechanisms predict that Trojans would survive 
the migration process, there are alternative models of planet migration that predict Trojans would 
not survive. For example, while planet formation models generally agree that planets should form 
on nearly circular orbits, it is possible that gravitational perturbations by other planets or a binary 
companion could excite sizable eccentricities. One possible formation mechanism for short-period 
giant planets is that a planet acquires a large eccentricity (e.g., due to strong planet-planet scat- 
tering, secular perturbations from a binary companion, or being tidal captured) and comes so close 
to the star that tidal dissipation circularizes the orbit at a small semimajor axis (Rasio &: Ford 
1996, Wu & Murray 2003, Gaudi 2003, Ford & Rasio 2006). Detecting a Trojan companion to 
a short-period planet would present a serious challenge for these mechanisms for forming "hot 
Jupiters" and would imply that the planet in such a system was formed via migration through a 
dissipative disk. Thus, searching for extrasolar Trojans can test models of planet formation. Here, 
we present a method for detecting Trojan companions to extrasolar planets by combining RV and 
photometric observations of transiting extrasolar planets. We refer to all bodies librating about the 
L4/L5 fixed point of a planet as "Trojans" and focus our attention on Trojans that are significantly 
less massive than the currently known planet and not currently recognizable from radial velocity 
(RV) observations alone. 
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2. Observational Constraints on Trojans 



We denote the stellar mass (m*), the planet mass (m p ), and the Trojan mass (my). Since the 
known transiting planets have short orbital periods and are subject to rapid eccentricity damping 
(Rasio et al. 1996), we initially assume the planet to be on a circular orbit about a star. Then a 
Trojan would orbit at one of the two fixed points, L4/L5, which lie along the orbit of the planet 
and lead/trail the planet by 60°. If there are no other bodies in the system, then the L4/L5 fixed 
points are stable if the ratio, ji = (m p + mT)/(m* + m p + iut), is less than a critical threshold /j, c , 
where 0.03812 < \i c < 0.03852 and \x c depends on e = rriT/(m p + mx) (Murray & Dermott 2000). 
If the Trojan resides exactly at the L4/L5 fixed point, then the direction of the vector sum of the 
forces exerted on the star by the planet and Trojan will lead/trail the force exerted on the star by 
the planet alone by an angle, 4>, such that tan</> ~ y/Se/(2 — e) x (1 + 0(/i)) (Fig. 1). 

More generally, for a Trojan that is librating about the L4/L5 fixed point, 4> will vary by an 
angle |A</>| ~ A^f ast + A s i ow > where A^f ast varies on the orbital period of the planet P, and A^ s i ow 
varies on the secular timescale, Put, ~ Pa/ 4/27 / u~ 1 / 2 (Murray & Dermott 2000). Since the Trojans 
of short period planets are likely to have formed in the presence of a dissipative disk, we focus on 
Trojans undergoing small librations with an amplitude, da <C / u 1 / 2 a p , where a p is the semimajor 
axis of the planet. In this case, A0f ast ~ eye, where ex is the Trojan's osculating eccentricity, 
and A</> s i ow ~ [S a/(6 / u 1 / 2 a)]e. Nesvorny et al. (2002) show that the behavior of Trojans of planets 
with small eccentricities is similar. Perturbations (e.g., GR, stellar quadrupole, tides) are unlikely 
to affect the stability, since they typically have timescales much longer than the orbital or libration 
timescale. 

If a planet on a circular orbit were the only body perturbing the central star, then the time 
that the stellar RV equals the RV of the system barycenter (To) would coincide with the time of 
midtransit (T c ). However, the gravitational perturbation of a Trojan at L4/L5 would cause these 
two times to differ by 



A Trojan can signal its presence by a time offset between the ephemeris determined from transit 
photometry and the ephemeris determined from RVs (Fig. 1). If there are Trojans at both L4 
and L5, then this measures the difference in mass at L4 and the mass at L5. If a planet is on 
a slightly eccentric orbit, then there is an offset of A t ~ P/(2ir) x (ecosu; + 0(e 2 )), where u> 
is the argument of pericenter, even in the absence of a Trojan. While short-period planets are 
expected to circularize rapidly, it is desirable to constrain the eccentricity observationally (e.g., 
by RV observations or timing of the secondary eclipse) before claiming the detection of a Trojan. 
Additional planets could also perturb the time of midtransit (Holman & Murray 2005; Agol et al. 
2005) so the offset will vary from transit to transit. Therefore multiple transits should be observed 
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to verify that any observed offsets are not due to perturbations by a more distant planet. 

For a transiting planet, both P and T c can be measured precisely using photometry alone. 
Consider a series of continuous photometric observations with uncorrelated Gaussian uncertainties 
of magnitude a p h, taken at a rate T around a single transit. The transit time can be measured 
with an accuracy of <tt c ~ ytj2Ta p hP~ 2 , where t e is the duration of ingress/egress and p is the 
ratio of the planet radius to stellar radius. For typical parameters (e.g., a p h ~ 10~ 3 ), T c can be 
measured to ~ 10s (e.g., Brown et al. 2001). The period can be measured much more accurately, 
from observations of multiple transits separated by many orbits. 

Given the measurement precision for T c and P, the practical limit on measuring At is set by 
the uncertainty in To from RV observations. Assuming a circular orbit, the RV observations at a 
time U can be fit by the mode0 v { = C + K sin (2n(ti - T c )/P + <j>) = C + A sin (2tt(^ - T c )/P) + 
B cos (2ir(ti — T c )/P). where A = K cos (ft, B = IT sin 0, and K is the velocity semi-amplitude. 
Assuming the period determined from photometric observations, the coefficients A, B, and C, (and 
hence the phase difference, tan</> = B/A) can be determined by linear least squares fitting to RV 
observations. If there are Nry RV observations with uncorrelated Gaussian uncertainties {(Try) 
and many RV observations are evenly distributed over orbital phase, then a Fisher information 
analysis (Gaudi & Winn 2006) reveals that the uncertainties in model parameters will approach 
a A = cr B = y/2/NaRv, = ^2/N rv or V /K, and a A t ^ o"T = V 1 / 2vr 2 N RV Pctrv / K . A similar 
analysis for an eccentric orbit in the epicyclic approximation, shows that the uncertainty in At 
is increased by a modest factor over the expression above. If a Trojan were present, then the 
uncertainty in <f> would set the uncertainty in the measurement of the mass of the Trojan, 
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If we were to demand a measurement of At > 3.291o"At to claim the detection of a Trojan, then 
a total of ~ 160 (60) precision RV measurements could detect a ~ 3M® Trojan, assum- 

ing a host star with an intrinsic jitter, Oj ~ 3 m s _1 (Wright 2005), and 1 m s~ 1 measurement 
uncertainties added in quadrature. While challenging, it is remarkable that current ground based 
instruments have the necessary precision to detect such a low mass Trojan with a plausible amount 
of observing time. 



ignoring any observations during primary transit when the Rossiter-McLaughlin effect distorts the observed RV 
(Winn et al. 2005; Gaudi & Winn 2006). 
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3. Example Application 

In Table 1, we summarize the current observational parameters and sensitivity to Trojan 
companions of extrasolar planets that transit bright stars, based on the above analysis. We find 
that combining the above method with existing observations already provides significant upper 
limits on the mass of Trojan companions to the planets HD 209458b and HD 149026b. Next, 
we perform more careful Bayesian analyses of the current observational constraints for these two 
cases. For HD 209458b we adopt the transit period and ephemeris of Knutson et al. (2006). We 
reanalyzed the RV measurements from Butler et al. (2006), fixing the orbital period and transit 
ephemeris. We use Markov chain Monte Carlo (Ford 2005, 2006; Gregory 2005) to sample from 
the posterior probability distribution for the remaining RV model parameters K, e, to, Mq, C, and 
Cj, where Mq is the mean anomaly at the epoch of midtransit. We assume priors that are fiat in 
log(l + K/K Q ), e, to, Mq, C, and log(l + aj/a Q ), and choose K Q = a Q = 1 m s _1 , but our results 
are insensitive to these assumptions. We then construct the posterior distribution for the quantity 
(Mo — ecosa;)P/(27r) ~ At. In Fig. 2a we show the distributions for At using three different 
assumptions. We find At = —11.4 ± 8.7 min (circular orbit), At = —16.4 ± 10.8 min (eccentric 
orbit ignoring the secondary eclipse), and At = —13.1 ± 8.9 min (eccentric orbit using the time of 
the secondary eclipse; Deming et al. 2005). We conclude that existing observations place an upper 
limit on the mass of Trojan companions to HD 209458b of 13.2M® at the 99.9% confidence level. 

We have performed a similar analysis of HD 149026b (Fig. 2b) using the observations of Butler 
et al. (2006) and Charbonneau et al. (2006). If we assume a circular (eccentric) orbit, then we find 
At = —19 ± 31 min (At = 98 ± 112 min). The constraint is significantly weaker when we allow 
for an eccentric orbit, due to the limited number of RV observations and poor phase coverage. 
Incorporating a preliminary estimate of the time of the secondary eclipse (J. Harrington 2006, 
private communication), we find A t = 13 ± 27 min and place an upper limit on the mass of Trojan 
companions to HD 149026b of 24.5M e at the 99.9% confidence level. 

4. Discussion 

In principle, Trojans could be detected via their radial velocity, astrometric, transit, or transit 
timing signatures. If a Trojan is sufficiently massive and has a sufficiently large libration amplitude, 
then it could be detected from the deviations from a Keplerian perturbation to the stellar radial 
velocity or astrometric signal caused by a single planet. Laughlin & Chambers (2002) have shown 
that two comparable mass planets occupying a 1:1 mean motion resonance would typically have 
strong planet-planet gravitational interactions on a secular timescale. However, these signatures 
may not be unique: a reanalysis of the RV observations of HD 128311 and HD 82943 have shown 
that both of the current data sets are consistent with a pair of planets in a 1:1 mean motion 
resonance (Gozdziewski &: Konacki 2006), as well as the originally published orbital solutions 
involving higher-order mean motion resonances. 
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Trojans may also be detectable if they transit their parent star. Photometric or spectroscopic 
monitoring of stars with transiting planets (particularly at times offset from the planet transit by 
~ P/6) may reveal the Trojan transit via the decrease in stellar flux or anomalous RV excursions 
due to the Rossiter-McLaughlin effect (Gaudi h Winn 2006). For planets discovered via a photo- 
metric transit search, there will typically be observations at epochs useful for searching for large 
Trojans. Unfortunately, a Trojan might not transit its parent star if it has a significant inclination 
(e.g., Morbidelli et al. 2005). Since the libration period can be large, long-term monitoring would 
be required to ensure detection. Currently, the most stringent photometric constraints on Trojan 
companions to HD 209458b come from the continuous photometry of the system for 14 days by the 
MOST satellite (Rowe et al. 2006). When heavily binned into ~ 2 hr intervals, this photometry 
has a fractional uncertainty of ~ 3 x 10~ 4 (Rowe et al. 2006). Assuming an average density equal 
to that of Earth, this corresponds to a 3-a detection threshold of 48M e . It is not clear whether the 
data reduction techniques used in their analysis of HD 209458 might subtract part of the signal due 
to a Trojan (Rowe et al. 2006). Regardless, if a Trojan had a vertical libration amplitude greater 
than ~ 9°, it would not always transit the star. Since the putative libration period of ~ 53 d 
(Murray &i Dermott 2000) is significantly longer than the duration of the MOST observations, it 
could have missed even a much larger Trojan. 

In a sense, our method is most similar and complementary to the recently-proposed method of 
searching for gravitational perturbations due to low mass planets using transit timing (Holman & 
Murray 2005; Agol et al. 2005). In contrast to the transit timing method, the unique geometry of 
Trojan orbits results in a nearly constant perturbation (assuming small amplitude libration about 
L4/L5) that recurs at every transit. In principle, it is not necessary to make precise measurements 
of the time of many transits to search for a complex pattern of perturbations. Thus, our method 
can be practically applied to transiting planets with long periods. Nevertheless, multiple transits 
should be observed to ensure that the same offset is observed and avoid potential confusion with 
perturbations from a more distant planet (Holman & Murray 2005; Agol et al. 2005). 

In principle, our technique could be applied to search for terrestrial-mass Trojans of giant 
planets orbiting in the habitable zone of their stars (Ji et al. 2005; Schwarz et al. 2005). While 
present search techniques are strongly biased towards finding transiting planets at short orbital 
periods, future space missions (e.g., Corot, Kepler) offer the prospect of finding transiting planets 
in the habitable zone of their stars, particularly for low mass stars where the habitable zone can 
be ~ 0.015 AU away from the star. 

We thank E. Agol, T. Beatty, E. Chiang, D. Fabrycky, M. Holman, R. Nelson, and J. Winn 
for helpful comments. We thank Joe Harrington for providing the time of secondary eclipse for HD 
149026b. Support for E.B.F. was provided by a Miller Research Fellowship and by NASA through 
Hubble Fellowship grant HST-HF-01195.01A awarded by the Space Telescope Science Institute, 
which is operated by the Association of Universities for Research in Astronomy, Inc., for NASA, 
under contract NAS 5-26555. Support for B.S.G. was provided by a Menzel Fellowship from Harvard 
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Table 1. Sensitivity to Trojans of Extrasolar Planets 



Star 


P 


K 


°RV & 


N RV 




<JAt b 


u iriT 


References 




(d) 


(m s- 1 ) 


(m s- 1 ) 




(M ) 


(min) 


(M e ) 




HD 209458 


3.5247455(2) 


84(1) 


5.0 


64 


1.101 


8.5 


2.6 


1,2 


HD 149026 


2.87598(1) 


43(2) 


5.7 


16 


1.30 


30.9 


6.2 


1,3,4,5 


HAT-P-1 


4.46529(9) 


60(2) 


5.1 


13 


1.11 


34.0 


6.4 


6 


TrES-2 


2.47063(1) 


181(3) 


6.9 


11 


1.08 


9.2 


7.6 


7 


HD 189733 


2.21857(2) 


205(6) 


15 


24 


0.82 


10.7 


8.9 


8,9 


TrES-1 


3.030065(8) 


115(6) 


14 


8 


0.89 


42.3 


17. 


1,10 


XO-1 


3.94153(3) 


116(9) 


15 


6 


1.0 


67.4 


25. 


11,12 



a When available, we list the rms velocity to the published best-fit RV model rather than the 
quoted measurement uncertainty. 

b We list "l-cr" uncertainties, implicitly assuming circular orbits for the transiting planets. 
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Fig. 1. — Illustration of the method to detect Trojan companions of transiting planets by comparing 
the transit and RV observations. Views of the star, planet and Trojan (A,C plan, B,D from 
observer's perspective; not to scale). The grey circle shows the orbit of the planet and Trojan. 
The dotted line indicates the direction of the acceleration of the star, the dashed line the direction 
of the transiting planet, and <j) 1S t ne angle between these two directions. The vector shows the 
direction toward the observer. Panels (A,B) show the position at To, the time of the stellar reflex 
RV null. Panels (C,D) show the position at T c , the time of the central transit. Panel (E) shows 
the stellar reflex RV as a function of time (in units of the period of the planet), with the times To 
and T c indicated. Panel (F) shows the intensity of the star as a function of time. We have assumed 
that the Trojan is inclined so that it does not transit the parent star. 
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Fig. 2. — Marginal posterior probability distributions of A t for HD 209458b (top) and HD 149026b 
(bottom). Here At ~ (Mo — ecosa;)P/(27r) is the difference between the time of the stellar reflex 
RV null and the time of central transit that could be due to a Trojan. The dotted curves assume 
a circular orbit, the dashed curves allow for a non-circular orbit (ignoring the constraint from the 
secondary eclipse) , and the solid curves allow for a non-circular orbit and incorporates the measured 
times of secondary eclipse. 



